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A model of two Calogero-Sutherland Bose gases A and B with strong odd-wave AB attractions 
induced by a p-wave AB Feshbach resonance is studied. The ground state wave function is found 
analytically by a Bose-Bose duality mapping, which permits to accurately determine static physical 
properties by a Monte Carlo method. The condensation of particles or particle pairs (molecules) is 
tested by analyzing the presence of the off-diagonal long-range order in one- or two- body density 
matrices. The p-wave symmetry of AB interaction makes possible quasi-condensation of type A 
particles at the Fermi momentum of the B component. The zero-temperature phase diagram is 
drawn in terms of densities and interaction strengths. 
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I. INTRODUCTION 

Strong interatomic interactions and correlations occur 
in ultracold gases confined in de Broglie waveguides with 
transverse trapping so tight that the atomic dynamics is 
essentially one-dimensional (ID) with confinement- 
induced resonances [l], allowing Feshbach resonance 
tuning 3] of the effective ID interactions to very large 
values. This has led to experimental verification [J-ly] 
of the fermionization of bosonic ultracold vapors in such 
geometries predicted by the Fermi-Bose (FB) mapping 
method Q, an exact mapping of a ID gas of bosons with 
point hard core repulsions, the "Tonks- Girardeau" (TG) 
gas, to an ideal spin-aligned Fermi gas. The "fermionic 
Tonks-Girardeau" (FTG) gas [1, @, a ID spin-aligned 
Fermi gas with very strong attractive interactions, can 
be realized by a 3D p-wave Feshbach resonance as, e_^ 
in ultracold "'"K vapor [T^]. It was shown recently 
that a mixture of two ID ideal Bose gases A and B with a 
strongly attractive AB interaction of the same FTG form 
is exactly solvable by a Bosc-Bosc duality mapping to a 
mixture of two ideal Bose gases with no AB interaction. 
Reduced density matrices of all orders were shown to 
be reducible to ID integrals, and it was found that the 
off-diagonal elements of the one- and two-particle den- 
sity matrices have unusual behavior: The strong AB at- 
traction destroys the ground state Bose- Einstein conden- 
sation (BEG) and single-particle off-diagonal long-range 
order (ODLRO) of both components A and B, and it 
induces both AA and BB pairing manifested in super- 
conductive ODLRO of the two-particle density matrices 
of components A and B, although there are no A A or BB 
interactions. Furthermore, there is no AB pair ODLRO 
in spite of the strong AB attractions. It was also shown 



that if the AB attraction is a finite odd-wave attraction 
rather than the infinite FTG limit, and there is also a 
repulsive even-wave AB interaction of Lieb-Liniger (LL) 
(5-function form [12[ , then there is a quantum phase tran- 
sition as the coupling constants are varied, between a 
phase in which there are no AB contact nodes and only 
the repulsive LL interaction acts, and another phase in 
which there are AB contact nodes and only the attractive 
FTG-like interaction acts. 

The Calogero-Sutherland (CS) model [13] of a ID 
Bose or Fermi gas with interparticle interaction poten- 
tial A(A — l)/a:fj is analytically solvable for all values 
of the coupling constant A. Here we will show that a 
mixture of two CS Bose gases A and B with an AB inter- 
action of FTG form is also analytically solvable by the 
same Bose-Bose duality mapping described in the previ- 
ous paragraph, which in this case maps the system to a 
mixture of two CS gases with no AB interaction. This 
has the great advantage that the single-component CS 
model reduces to previously known results in several lim- 
its, thus providing an exact interpolation between those 
limits: ideal Bose gas for A = 0, TG gas for A = 1, and 
the limit of a classical crystal A oo. For the mixture, 
the ground state can be found in closed form even if com- 
ponents A and B have different CS coupling constants Aa 
and Ab. The ground state will be found in explicit an- 
alytical form and the one- and two-particle density ma- 
trices will be determined numerically by a Monte Carlo 
method. 



II. SINGLE-COMPONENT 
CALOGERO-SUTHERLAND MODEL 
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Before studying the properties of mixtures it is useful 
to recall some properties of a single-component CS sys- 
tem. For N particles of mass to in a periodic box of size 



2 



L, the CS interaction potential is 



TT^h^ A(A-l) 



(1) 



where A > is the interaction parameter. This interac- 
tion potential satisfies periodic boundary conditions and 
can be obtained by evaluating the L-periodic extension of 
an inverse square potential: V'~"^{x) — J2'jL-oa f^^H^ ~ 
l)/[m{x + jL^]. Another physical interpretation of the 
sine function in Eq. ([!]) is that particles stay on a ring of 
diameter L but the two-particle interaction corresponds 
to the chord distance Lsm{TTx/L). In the thermody- 
namic limit the interaction potential reduces to an in- 
verse square potential. Such a potential is very special 
as it scales as the kinetic energy and there is no any other 
length scale in the system different from density. Prop- 
erties of the system are then governed by the interaction 
parameter A. The bosonic ground-state wave function 
of particles interacting with the potential Eq. ([T|) can be 
written explicitly: 



N 

i;'^^{xi,...,XN)=Y[ 

i<3 



Tr{xi 



L 



(2) 



An important observation is that the tails of the inter- 
action potentials vanish for two special values A = and 
A = 1. In the first case the ground state wave function 
reduces to a constant which is the solution for an ideal 
Bose gas, while in the second case the solution ([2]) coin- 
cides with the absolute value of the wave function of an 
ideal Fermi gas and corresponds to the TG gas 

Although at present there are no realizations of the 
CS system, this model describes many important physi- 
cal regimes that can be reached in other one-dimensional 
systems. Furthermore, the long-range properties are ex- 
pected to be similar. In order to see that we note that 
the long-range wave function of a one-dimensional system 
can be deduced from zero-point motion of the phonons. 
The ground state has the form [l^ 'ip{xi, ...,xn) = 
V'short(a;i, ...,a;7v)V'phon(a;i, ...,a;7v) where V'shoit describes 
the short-range behavior and 
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Ipphonixi, ...,Xn) Y\. 

i<j 



TT{Xi -Xj) 



L 



i/2 



(3) 



for large separations between particles. Here a is a con- 
stant related to the speed of sound. Comparing ^ to 
([2]) we see that actually the CS wave function retains the 
structure typical of phonons even at short distances. 

The presence of phonons in a one-dimensional sys- 
tem permits construction of an effective Luttinger liquid 
model, which is based on the assumption that the exci- 
tation spectrum is linear in momentum Ek = h\k\/2mc, 
where c is the speed of sound. The Luttinger liquid yields 
the long-range behavior of the correlation functions and 
leads to some important conclusions on the generic be- 
havior of one-dimensional systems, such as the absence of 



Bose-Einstein condensation and absence of crystalline or- 
der. Within this approach properties are defined by one 
governing parameter, the Luttinger constant K — vp/'ic 
with vp — irhn/m being the Fermi velocity. The speed of 
sound in the CS model can be calculated from the com- 
pressibility mc^ — ndfi/dn (as usual fi = dE/dN is the 
chemical potential) and leads to a very simple relation 
K = 1/A. This means, in particular, that the CS model 
has the same Luttinger parameter as the (5-interacting LL 
gas in the region < A < 1. Values of A slightly larger 
than unity correspond to the "super- Tonks-Girardeau" 
regime. Finally, A > 2 correspond to a quasi-crystalline 
region, with diverging peaks in the static structure fac- 
tor. From this point of view the CS model describes 
long-range properties in very different physical regimes, 
ranging from the ideal Bose gas to a classical crystal. 
For a detailed study of these various regimes see [l6| . 
In the regime A < 1/2, in addition to the well-behaved 
solution ^ the CS model possesses singular solutions as- 
sociated with "fall to the center" [13, El]- Nevertheless, 
the state ^ remains well-behaved; we conjecture that 
for < A < 1/2 it is the lowest gas-like state. The long- 
range properties of this state agree with those of the LL 
gas in the whole regime < A < 1. 



III. FTG INTERACTION 

The FTG gas is a spin-aligned ID Fermi gas with in- 
finitely strongly attractive zero-range odd-wave interac- 
tion induced by a p- wave Feshbach resonance. It is the in- 
finite ID scattering length limit aiD — oo of a ID Fermi 
gas with zero-range attractive interactions leading to a 
ID scattering length defined in terms of the ratio of the 
derivative of the wave function to its value at contact: 
^'(xjfc = 0+) = -'i'ixjk = 0-) = -aiu'i'ixjk = 0±) 
where the prime denotes the derivative with respect to 
Xjk @, BQ- The FTG limit aiD —oo corresponds to 
a ID zero-energy odd-wave scattering resonance reach- 
able by Feshbach resonance tuning to a ID odd-wave 
confinement-induced resonance [ij-yl . An AB interaction 
of this form was used in Ref. |11| for a mixture of two 
ID ideal Bose gases A and B, and here we will use it for 
the AB interaction between two CS Bose gases A and B; 
Ref. [ll| can be consulted for more details of the inter- 
action and its effect. 



IV. TWO-COMPONENT MIXTURE AND 
MAPPING SOLUTION 



Consider now a mixture of two CS gases A and B with 
an AB interaction of FTG form and no longitudinal trap 
potential, with periodic boundary conditions of periodic- 
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ity length L. The Hamiltonian is 



^ 2mA dx^ 

1 = 1 * i<i<J<NA 

N-B 



^=l j=l 



(4) 



where (a;i, • • • ,xnj^) are the particle coordinates for com- 
ponent A, {ui,--- tUn-b) ^-re those for component B, 
Va{x) and Vb{x) are the previously-defined CS interac- 
tion V'^^{x) with CS coupling constants Aa and Ab, and 
^'FTG (x) is the FTG interaction described in 0, H, S IHl 
and Sec. Ill above. The exact ground state of the model 
of [ll[ was found with the aid of a Bose-Bose duality 
mapping function [l^ 

Na Nb 

M{xi,--- ,XNA;yi,--- ,yNB) = YlYisgnix.-yj) , (5) 

1=1 J=l 

where the sign function sgn(a::) is -1-1 (—1) if a; > (a; < 
0), and the same mapping can be used to obtain the 
exact solution of the present model. Let be the exact 
ground state of our mixture of two CS gases A and B with 
FTG AB interactions, which we wish to determine. As 
in flT], mapping by AI removes the interaction between 
components A and B, in the sense that = '^mqM 



where Vfivio, the ground state of a mixture of two CS gases 
A and B with no AB interactions, is the product ^"1^10 — 
^OA^OB of the ground states of CS gases A and B. The 
explicit expression for the ground state wave function is 

Na Na Nb Nb 
i=l jyi k=l l>k 

I sin7r(yfc - yij/Ll^"" sgn{x, - yu) ■ (6) 

Since [M{xi,--- , xat^ ; yi, • ■ ' ^UNb)? = 1: all diagonal 
density matrix elements of the ground state '^0 with 
FTG AB interactions are the same as those of the model 
ground state which has no AB interactions. In 

particular, the A and B component densities in the 
present case (no trap potential) are trivial constants, 
uaIx) = Na/L and nB{y) = N-q/L where L is the pe- 
riod of the periodic boundary conditions. Furthermore, 
the AA and BB pair distribution functions are just those 
of single-component CS gases with coupling constants Aa 
and Ab, and the AB pair distribution function is just a 
constant nA«B- 



V. ONE-PARTICLE DENSITY MATRICES AND 
MOMENTUM DISTRIBUTIONS 

Next consider the off-diagonal elements of the single- 
particle density matrices Pia{x, x') and pisiy, y') of com- 
ponents A and B. One has 



PIa{x,x') = NaL ^* J 'i>o{x,X2, ■ ■ ■ ,XNA;y)'^o{x',X2, ■ ■ ■ ,XNAiY)dX2 ■ ■■dxNAdY 
n Nb 

= / ^'oA(a;,a;2,-- - ,xna)'^ok{x' ,X2, - ■ ■ , a:^7VA)*0B(yir ■ • ,yNB)\\_^S^{x - yj)sgn{x' - yj)dx2 ■ ■ ■ dxNAdY 

„ Nb 

\x,x') J 4'ob(2/i>-- - , VNb) Y[sgn{x - yj) sgn{x' - yj)dY , (7) 



PIT' 



where Y = (j/i, • ■ • , j/atj,) and p^™^{x, x') is the single-particle density matrix of the pure (single-component) CS gas. 
The formula for piB(y,y') is the same, with obvious interchanges of A with B and x with y. 
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VI. TWO-PARTICLE DENSITY MATRICES 



By generalization of ([7]) and the mapping used in [ill [20| the two-particle A-component density matrix 
P2AA(a;i,x2;xi,x2) = Nk{Nk - 1)L~^'^^+^^^ J ^o{xi,x2,X3,- ■ ■ , xtv*; ^)*o(2;'i, a;2, 2:3, • ■ • , xata; i^)da;3 • ■ • rfa^WA^^^ 

*OA(a;i,--- ,XNf,)'if0A{x'i,x'2,X3, - ■ ■ ,XNA)'^OBiyir ■ ■ ,2/7Vb) 

]J sgn(a;i - yj) sgn{x[ ~ yj) sgn(a;2 - y^) sgn(a;2 - yj)dx3 ■ ■ ■ dxN^dY 



Nb 



„ Nb 

P2Aa(^1' 2^2; x'l, x^) / *^B(yi, ■ • • , yjvB) H '^gn(xi - y^) sgn(x'i - y^) sgn(a;2 - y^) sgn(4 - yj)dY 



(8) 



where /52Aa(^1i ■^2! a;^, 0:2) is the two-particle density matrix of the pure CS gas, and the formula for p2BB differs only 
by the obvious interchanges. The AB-pair density matrix is 



P2AB(a;,y;a;',y') = NaNbL (^a+JVb) y 'q)Q[x,X2,X3, ■ ■ ■ , a^Ar^; y, y2, • • • , yAfB)*o(a;', 2:2, • • • ,XNj^;y' ,y2, - ■ ■ ,VNb) 

X dx2 ■ ■■dxNf,dy2 ■ • • dyjVB = j '^0a{x,X2, ■ ■ ■ ,XNA)'^0B{y,y2, ■ ■ ■ ,yNB)'^OA{x' ,X2, ■ ■ ■ ,XNA)'^0B{y',y2, ■ ■ ■ ,yNB) 

X sgn(a:: - y) sgn(a::' - y') 



n ^g^(^ ~ yj") - Vj) 

J=2 



■iVA 



Yl sgn{x^ - y) sgn(a::i - y') 



.i=2 



dx2 ■ ■ ■dxNAdy2 ■ ■ -dyNB 



(9) 



VII. MONTE CARLO TECHNIQUE 

In order to obtain numerically the correlation func- 
tions we resort to Monte Carlo methods. Such meth- 
ods are very efficient for evaluation of multidimen- 
sional integrals which in our case correspond to op- 
erator averages {A) / ... / A(a;i, a^A^A, yi, y^B) 
|4'o(a:i,.^..,a;ArA,yi,...,yArB)|^ dxi ... dx^^ dyi ... dyAfB, 
where A is some operator of interest. We obtained 
an explicit expression for the ground state wave func- 
tion vl'o, as given by Eq. (jG)), which significantly simpli- 
fies the numerical calculations. The basic idea of the 
method is to generate a sequence of points Ri in the 
phase space R — {xi, ...,X]\jj^,yi, ...,yN-g} according to 
the probability distribution |5'o(R)p. Then the opera- 
tor average is approximated by an average over R^ as 
(A) « J2i=N^,asur-. M'^t) /^measure- The Statistical er- 
ror of the estimation is kept under control and can be 
reduced by increasing the series of measurements. We 
follow Metropolis prescription '21'] for generation of the 
sequence of R^ such that each new point in the phase 
depends only on the previous point Ri_i (i.e. such a 
sequence forms a Markov chain): 



• generate a new configuration R.; 

• the move from Ri_i to Ri is always accepted if the 
wave function in the trial configuration is larger, 

*o(R») > *o(R»-i) 

• if the wave function in the trial configuration 
is smaller, accept the move with the probability 
[vI'o(R.)/*o(R.-iF 

We generate new moves by displacing a random particle 
by a distance taken from a Gaussian distribution. The 
width of the Gaussian is adjusted in such a way that the 
acceptance rate is close to 50%, apart from the singu- 
lar case of vanishing Aa and Ab, when the probability 
distribution degenerates |\l/(R)p — 1 and any move is 
accepted. In this case we adjust the width of the Gaus- 
sian in such a way that the typical displacement is small 
compared to the system size. 

The one-particle density matrix ([7]) can be easily re- 
cast in a form suitable for Monte Carlo calculation. In- 
deed, taking into account that in a homogeneous system 
Pia{x, x') is a function of the difference of arguments and 
labeling this difference as x one has 
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riA 



J dXdY 



(10) 



The interpretation of formula (jlOp is that the The physical importance of p2AA(a;i, a;2; x'j^, is that it 
ratio of values of the wave function [5'o(a;i + 
x,X2,--- ,XNj^;Y)/"i'Q{X;Y)] with and without a par- 
ticle displaced by x is averaged over the random walk 
distributed according to ^q(X, y). 

Most generally, the two-body particle density ma- 
trix ([5]) depends on four variables, namely, xi,X2,x[,X2- 

I 



can be used for testing a possible formation of a Bose con- 
densate of particle pairs (molecules). One should study 
its behavior when a pair xi, X2 is displaced by some dis- 
tance R, i.e. x'l = xi + R and x'2 = X2 -\- R- Then the 
relevant parameters are the displacement length R and 
the size of the pair r ~ xi — X2 — x'l ~ x'2. 



P2AAir, R) 



J_\ ![Six2 -xi~ r)^o(a;i + R,X2 + R,X3,- ■ ■ , xn^;Y)/^q{X; Y)]^liX; Y)dXdY 
Na) j^l{X-Y)dXdY 



(11) 



r 



The ODLRO of pairs is manifested by non-zero large- 
R asymptotic value of the projected TBDM defined 
from ([TT|) by integrating out the pair size r|22|-|24| 



P2{R) 



2 



dxidx2P2{xi+ R.,X2 + R,xi,X2) ■ (12) 



If ODLRO of pairs (molecules) is present, the two-body 
density matrix at large separation distances R reduces to 
the product of molecule orbitals. In Monte Carlo calcu- 
lations this is tested by setting i? to L/4 in Eq. pT|) . 



VIII. RESULTS 

A. Mixture of ideal Bose gases with FTG 
interaction 

We will consider first the case of vanishing Aa = Ab = 
0. This limit describes a system consisting of a mixture 
of two ideal Bose gases with FTG interaction between 
different components. Such a system is interesting as for 
a single component CS system A = is the only value 
of the interaction parameter for which true Bose conden- 
sation exists, see Section [TTl Instead for any finite A the 
OBDM has zero asymptotic value. As well it is conve- 
nient to start the discussion from the case of mixtures of 
ideal Bose gases, as this the limit where explicit analyt- 
ical expressions for the one- and two-particle correlation 
functions are known from Ref. |1J| . 

In the absence of the other component, ?ia = 0, the 
OBDM is constant Pia{x) = ha and the Bose condensa- 
tion is complete. For a finite concentration of the second 
component, the one-body density matrix can be obtained 
explicitly even in a system of a finite size and is given by 
Pia{x) = nA(l — 2nB|a;|/7VB)^^. In the thermodynamic 



limit the decay in pia{x) is exponentially fast 

Pia(x) = nAe-2"-l^l (13) 

and the strength of the decay depends on the density of 
the other component ub- In particular, for zero concen- 
tration the strength of the decay is zero and the result 
for the ideal Bose gas is recovered. Instead, any finite 
concentration of the other component removes the Bose 
condensate in a strong sense, that is exponentially fast 
decay compared to power-law decay in interacting one- 
dimensional Bose gas (see Fig. [T]). 

The momentum distribution is related to the one-body 
density matrix (jl3p by a Fourier transformation and has 
a Lorentzian shape [il[ 



"AlKj - .2 , 



4ng -t- fc2 



(14) 



An important feature of n{k) in this case is that the mo- 
mentum distribution is finite for zero momentum with 
its value defined by the ratio of the densities nA(0) = 
n-A/nB- Another important observation is that the high- 
momentum decay follows the law 1/fc^ contrary to l/fc"* 
decay of a Lieb-Liniger of Tonks-Girardeau gas. We will 
comment more on the high-momentum behavior later. 

It was shown in Ref. ^] that while FTG interaction 
destroys Bose condensation in each component, at the 
same time it induces ODLRO in TBDMs p2AA and p2BB • 
Indeed, it was found that the spectral representation of 
the TBDM 

P2ixi,X2;x[,x[) = ^ Xi(t)*{xi - X2)(l)i{x[ - X2) (15) 

i 

reduces in the physically interesting regime Xi,X2 ^ 
x'i,X2 (for example, this is the case when a pair is dis- 



6 



1.0 



f 0.8 



0.6 



0.4 



0.2 



0.0 



-0.2 



exp(-2n|,|x|) 




>.=0 


?.=0.6 


X=0.1 


1=0.7 


^=0.2 


X=O.S 


?.=0.3 


>.=0.9 


X=OA 


l=^ 


?.=0.5 


x=^o 



100 



1 2 3 4 5 

FIG. 1: One-body density matrix pia{x) in a balanced system 
TVa ~ Nb ~ 10 for different interaction strength Aa = Ab- 
Thin uppermost curve, Aa = Ab = (mixture of two ideal 
Bose gases), thermodynamic limit as given by Eg. 1131 Thick 
lines, Monte Carlo results for A'a = A^b = 10 particles and 
increasing interaction strength (from upper to lower curve) 
Aa = Ab = 0; 0.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9; 1; 10. 



placed by a large distance) to 

P2Aa(^i,X2;x'i,4) = nie-2"Bki-x2|g-2„B|x;-4l . 

(16) 

In this way the largest eigenvalue A — Np^nj^/2n-B is 
macroscopic and the corresponding eigenfunction 



■2nB |a; I 



(17) 



can be interpreted as a dimer orbital. There is a BEC- 
BCS-like crossover from AA-pair EEC when n-g, ^ ua 
and the range of (j)^i is <C I/ua, implying tightly 
bound AA pairs, to extended and strongly overlapping 
AA Cooper pairs when ne ^ ua and the range of 0^"^ 
is > l/uA- The TBDM P2Bb(x) of the other com- 
ponent exhibits similar behavior with A and B inter- 
changed. It is interesting to note that largest eigenvalue 
A = A^a^-a/S^b in a balanced Na ~ Nb = N system 
reduces to N/2, i.e. all N/2 dimers are condensed. In an 
unbalanced system the range of a wave function of AA 
dimer (|17p is governed by the density hb of the other 
component. The eigenfunctions have to be normalized 
to unity |(/>(a;)|^ Ldx = 1 and as a result the normal- 
ization constant of (f>^'^{x) depends on the density ub- 
In this way the ratio ha/tib appears in the eigenvalue 
rescaling Na/2 to larger or smaller values, according to 
the considered densities. 

In terms of the projected TBDM (i?) (see Eq. p^ ). 
the presence of ODLRO means a finite asymptotic — > 
oo value. Contrary to dimensionless OBDM, where 
ODLRO is manifested as a constant independent of the 
number of particles, in dimensionless TBDM the finite 
value decreases with the number of particles as 1 /N. One 
way to understand this is that in pi{x,x') a particle is 
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FIG. 2: Inset: projected two-body density matrix P2Aa(-R)i 
Eq. (|12|l , in a mixture of two ideal Bose gases for A'a = Nb = 
100; 80; 60; 40; 20; 10; 5 (from upper to lower curves). Main 
figure: the same function multiplied by 2N. 
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FIG. 3: Two-body density matrix p2AA{r, L/A) in a balanced 
system of A'a = Nb ~ 10 particles for different interaction 
strengths Aa = Ab = 0; 0.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9; 1 
(decreasing from the value at zero). In the presence of BEG 
of molecules p2AA{r, L /4) is proportional to the square of 
the molecular orbital. For A = the thermodynamic result 
[i't^iir)]^ iflT)) is shown for comparison. 



annihilated at x and created at x' . The value of the wave 
function remains the same if this particle belongs to the 
condensate. In p2{xi,X2\x'i,x'2) a pair is displaced. In 
order to find a constant for a given particle xi from the 
molecular condensate the particle X2 should belong to 
the same molecule. The probability of that is 1/{N — 1). 
We account for this by multiplying p|'(i?) by 2N in the 
main figure [2l and a constant value is observed. The nu- 
merical error increases with the number of particles be- 
cause the error is multiplied by 27V, and in addition the 
phase space becomes larger and Monte Carlo simulation 
becomes more time consuming. 

We have checked the dependence of the TBDM 
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FIG. 4: Schematic phase diagram for Aa = Ab = showing 
a crossover as the density imbalance na/nA is changed from 
zero to infinity. Component A passes through three phases 1: 
true BEC of A, 2: quasi-condensate of A, BEC of large A A 
dimers, 3: BEC of highly localized dimers (molecules). Com- 
ponent B passes through the same phases but in an inverted 
order. The system is normal in all regimes as no sound exits. 



P2{R,r) on R, i.e. on the distance to which a pair is 
displaced. Its finite asymptotic value for — > oo mani- 
fests the presence of ODLRO of molecules. For \R\ large 
enough, so that such an asymptotic value is reached, the 
shape of p2{R,r) in the r direction corresponds to the 
square of the wave function of the bound state (molecu- 
lar orbital) [(t>^i{r)]'^ ■ Fig. [3] shows the dependence on 
r in a system of A'a = Nb = 10 particles for a value 
R = L/A within the asymptotic regime. (For the depen- 
dence on R refer to Fig. [51) Importantly, the range of 
the molecular orbital of an AA pair is governed by the 
density tib of the other component. 

The schematic phase diagram is drawn in Fig. |4] in 
terms of the density imbalance tib/tia- One way to fol- 
low it is to fix the number of A particles and change the 
number of particles of the other component. When B 
particles are absent the model reduces to an ideal Bose 
gas of A particles. Here a true Bose condensation of 
A particles is predicted. The transition to the phase 
with finite riB is not continuous which makes the case 
7iB — an excluded point. Indeed, for any finite concen- 
tration of the other component, the asymptotic value of 
the OBDM is strictly zero. Still, the decay could be 
extremely slow and the phase coherence might be pre- 
served at distances large compared to the AA interparti- 
cle distance. We call this regime a "quasi-condensate" of 
A particles. Note that this notation differs from a qua- 
sicondensation in one-component systems (for example, 
Luttinger liquids), which in such systems is commonly 
understood as a slow power-law decay of OBDM. A fi- 
nite density of the other component riB > induces con- 
densation of A A dimers. The dimers are shallow in the 
regime of quasicondensation of A. Instead for tib 3> tia 
such dimers have extremely localized wave function (or- 
bital) and we interpret this regime as a Bose-Einstein 
condensation of AA molecules. 

The mapping preserves the excitation spectrum. From 
this we infer that the gas with Aa = As = is not 
superfluid, as no sound exists and the system is not com- 
pressible. 



B. A mixture of Calogero-Sutherland gases with 
interspecies interaction: Fermi-momentum Bose 
condensate 



In this Section we discuss properties of a two- 
component gas with Calogero-Sutherland intraspecies 
potential and FTG interspecies interaction. The ground 
state wave function is given explicitly in Eq. Al- 
though the Ideal Bose gas case is recovered in the limit 
-^A = -^B = 0, properties of the system at a small but 
finite value of A are qualitatively different. Indeed, the 
system becomes compressible due to interactions in each 
component. Appearance of the sound mode leads to su- 
perfluid behavior in a weakly interacting regime which 
has to be contrasted to a normal behavior of an ideal 
Bose gas. 

The p-wave Feshbach resonance which leads to the 
FTG interaction causes exotic properties which are quite 
uncommon for a bosonic system. Indeed, it seems natural 
to expect that the ground-state wave function of a con- 
ventional bosonic system is positive-definite and it has 
no nodes. But it is not the case for our system. Indeed, 
the ground state wave function ^ is symmetric under 
exchange of two A-A bosons or B-B bosons, according to 
bosonic statistics. As particles A and B are distinguish- 
able, no symmetry or antisymmetry is required with the 
respect to A-B exchange. This means that both s-wave 
and p-wave scattering channels are available. Commonly 
s-wave scattering is considered, although in a vicinity of 
a p-wave resonance the relevant scattering corresponds to 
a locally antisymmetric solution. As a result the many- 
body wave function vanishes when two particles meet and 
also it changes sign. 

The "Feynman's no-node theorem" !28i| applies to al- 
most all physical realizations of bosonic systems. The 
standard reasoning is that if the ground state solu- 
tion ^'o(xi, xat) has some nodes, its absolute value 
|\E'o(a;i, xn)\ is also a solution. Smoothing the kinks of 
the solution close to the node leads to lowering the kinetic 
energy, while the potential energy increases only slightly. 
From this one concludes that the ground state wave func- 
tion should be positive-definite and non-generate. This 
reasoning is not applicable in our case as the nodes in the 
wave function are due to the AB interaction potential and 
can not be removed. As a result the zero-temperature 
wave function is not positive definite and there are other 
solutions for the Hamiltonian that give the same ground 
state energy. 

The possibility of bypassing Feynman's argument is 
already present in the literature. The proposals are to 
use the long-lived metastable states of bosons in the high 
orbital bands of optical lattices as a result of "orbital- 
Hund's rule" ^29] and multi-component bosons wit h sp in- 
orbit coupling linearly dependent on momentum [30|. It 
is noted that the emergent states might experience un- 
conventional Bose-Einstein condensation at non-zero mo- 
mentum [3l| . It is proposed that a mixture of two-species 
bosonic atoms interacting through a p-wave Feshbach 
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resonance exhibits a finite-momentum atomic-molecular 
superfluid}32i]. 

Negative regions of the wave function will not change 
the properties of the local quantities (for example, AA 
pair distribution is the same as in a single component 
Calogero-Sutherland model), but dramatically modify 
non-local properties (for example, the one-body density 
matrix and the momentum distribution). In a single com- 
ponent gas with interaction parameter A the one-body 
density matrix has a power-law decay at large distances 
Pi{x) oc l/|ria;|'^/^ (all terms of the series expansions are 
explicitly obtained in Ref. [l^ using replica method) . For 
any finite value of A the asymptotic value of pi{x) van- 
ishes and the ODLRO in one-particle density matrix is 
absent. Still the concept of a quasi-condensate can be 
applied in the regime of small A, where the character- 
istic spreading of the one-body density matrix is large 
compared to the interparticle distance. 

In the one-dimensional system some analogies between 
fermions and bosons with p-wave interaction can be 
drawn. By fixing positions of all particles but one we 
expect to see oscillations in the sign of the wave func- 
tion as the particle is moved. In a Fermi system the 
sign changes each time the reference particle crosses an- 
other particle. In our system the sign changes each time 
A particle crosses B particle. The number of sign al- 
terations is proportional to the density of the particles 
that produce change of the sign. In a fermionic system 
the characteristic period of space oscillations is fixed by 
the Fermi momentum, which itself is proportional to the 
density kp = 7m. In our system the period of oscillations 
in one-body density matrix of A particles is fixed by the 
density of B particles with the corresponding momentum 
kp = ttub- We see that this characteristic momentum 
coincides with the Fermi momentum that the system B 
would have if it were fermionic. 

The momentum distribution nA{k) is shown for some 
characteristic values of A in Fig. [S] For a mixture of 
two ideal Bose gases the momentum distribution is a 
Lorentzian (jl4p with the maximum at zero momentum. 
In this case the value at the peak is equal to the den- 
sity of the corresponding component, i.e. nx{k) — ua- 
For A = 1 the maximum is shifted to the Fermi mo- 
mentum of the other component k — kp. For very 
large values of A there is again only one peak situated 
at zero momentum. In this regime the potential is much 
larger then the kinetic energy which leads to creation of 
a quasicrystalfl^. The momentum distribution can be 
approximated by a Gaussian, similarly to classical sys- 
tems. The Gaussian width can be extracted from the 
short-range series expansion of the one-body density ma- 
trix giA{x) = l—C2{nAx)'^ + ... by approximating it with a 
Gaussian giAix) = exp{— C2(»T.Aa;)^} j25] and calculating 
its Fourier transform 

^(fc)= /_cxp ^ . (18) 

The coefficient C2 is related to the kinetic energy, which 




FIG. 5: Momentum distribution for a balanced system ua ~ 
UB for characteristic values of the interaction parameter A — 
0; 1; 10. The solid line, A = 0, shows the thermodynamic 
limit of a mixture of two ideal Bose gases (|14p . Dashed line, 
A = 1, and dash-dotted line, A = 10 are results of Monte 
Carlo calculation with = 40 particles. Thin solid line, 
Gaussian approximation to the momentum distribution in the 
"classical" regime as given by Eq. (|18|) . Note that the height 
of the peak for A = 1 depends explicitly on the number of 
particles in the system. 




FIG. 6: 3D plot of the momentum distribution for a balanced 
system ua = wb as a function of the momentum and the 
interaction parameter A. 



itself can be obtained by using the Hellman-Feynman the- 
orem from the total energy. The resulting expression is 
C2 = AV(6(2A-1)) ;i^. 

The large-A: asymptotic of the momentum distribution 
is governed by the short-range physics and it can be ob- 
tained from the short-range behavior of the two-body 
Bijl-Jastrow terms. The Calogero-Sutherland interaction 
introduces a cusp when two particles meet each other as 
described by \sm{TT{xi — Xj))/L\^ term in the ground 
state wave function. The momentum distribution can be 
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written in the following form: 

n{k) — n j \^Q{ki,X2, ■■■)\'^dx2---dxN , (19) 

where we use the momentum representation of the 
ground state wave function with respect to the first 
particle (33} 

rL/2 

^oiki,-)= e-'^'^'^o{xi,...)dxi . (20) 

J-L/2 

Here for simplicity we write down only the first argument, 

i.e. ^'o(fci,...) = *o(fci,a;2,---,a;ArA;2/i;---:yA'B)- First 
we consider the case of the one-component Calogero- 
Sutherland gas, A^b = 0. For small distances between 
two particles the sine function in can be expanded 
to Tr\{xi — Xj)/L\^ with the relevant part of the integral 
being proportional to 

^o{k) oc J e-''''=\x\^dx . (21) 

As the relevant part of the integral is coming from small 
X we are free to change the limits of the integration and 
make them infinite. In order to make the resulting in- 
tegral convergent we multiply the integrant by an expo- 
nent exp(— a|a;|), evaluate the integral and then take the 
limit of a — )■ 0. This results in a power-law dependence 
^'o(fc, ...) cx The momentum distribution ([19]) 

then follows a power-law n{k) oc l/jfcp^-'^"'"'*'^ with the ex- 
ponent depending on the value of A. In a two-component 
system, A^b > 0, the integrand in (^0)) changes its sign 
each time x passes through a particle of the other type y^. 
This splits the integral (PT|) into different regions of inte- 
gration with the leading contribution from the each point 
Ui being of the order of l/|fc|. This changes the asymp- 
totic behavior of the momentum distribution behavior to 
n{k) (X 1/k'^ for all physical values of A > 0. As already 
noted above, the momentum distribution in the A = 
case can be explicitly evaluated as given by Eq. (IT4l) and 
it has a 1/fc^ tail. The case of A = 1 corresponds to 
the TG interaction potential with the Bijl-Jastrow term 
proportional to the absolute value of the interparticle 
distance \xi — Xj\. Similarly the Lieb-Liniger gas has 
a two-particle solution of \xi — xj — aiol, where oid is 
the one-dimensional s-wave scattering length. From this 
we conclude that the asymptotic of the momentum dis- 
tribution of TG and Lieb-Linger gases are n{k) oc l/|fc|'' 
for a single component system and n{k) oc l/|fcp for a 
two-component system with FTG interspecies potential. 

An example of finite-size dependence of the momen- 
tum distribution is shown in Fig. [T] The behavior of 
the tail is best seen on a logarithmic scale which we 
adopt in the main plot. The analytic prediction l//c^ 
is shown with a dashed line which on the logarithmic 
scale is parallel to the numerical curve, while the constant 
shift is related to the proportionality constant Casympt in 

n{k) = Casympt / k ■ 




FIG. 7: Inset: Momentum distribution of a balanced sys- 
tem with Aa = Ab = 1 for diflerent system sizes A^a = 
TVb = 10; 20; 40; 80; 160; 320; 1000 as a function of momen- 
tum in units of Fermi momentum k/kp. Main figure: the 
same data points sliown on a log-log scale as a function of 
the separation from the divergency point \k — fcpl/fcp. High- 
momentum l/k^ and low-momentum l/|fc|^''^ asymptotics are 
shown for comparison. 



The Fermi-momentum condensation peak is the most 
pronounced for values of A = 1 as can be seen from 
Fig. [6l We report the finite-size dependence of momen- 
tum distribution exactly for this value of A in Fig. [71 The 
one-dimensional analogue of Bose-Einstein condensation 
at some value of momentum k is understood as diver- 
gency of the momentum distribution at this momentum 
n(k) — >■ (X) in the thermodynamic limit N 00. The 
momentum distribution for difi^erent number of particles 
is shown on a linear scale in the instet of Fig. [Hi One sees 
that there are high peaks at values of momenta k = ±k-p 
and the height of the peaks increases with the number 
of particles. In order to find the functional form of the 
divergency it is convenient to plot the momentum distri- 
bution as a function of \k—kp \ on a log- log scale, as shown 
in the main plot of Fig. [6] We see that the divergency 
is well described by a power-law n{k) oc l/|fc|^/'* and the 
height of the peak increases as N'^/^. This should be con- 
trasted with the infrared divergency of a one-component 
Tonks-Girardeau gas of N^^^. 

For the Aa = 1 case which leads to the highest peak 
in nA(fc) at the Fermi momentum, we study the depen- 
dence on the BB interaction strength Ab . The results are 
shown in Fig. [8| When the other component is an ideal 
Bose gas, Ab = 0, the positions of B particles are uncor- 
related and the one-body density matrix giA{x) decays 
in a monotonic way to zero. In this case the peak in the 
Fermi momentum is absent. The oscillations in giA{x) 
are enhanced by ordering in the B component, which ap- 
pears in the system as Ab is increased. In particular, 
the height of the peak in nx{k) at the Fermi momentum 
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FIG. 8: Momentum distribution of a system with A'a = A^b = 
10 particles with the AA interaction strength fixed to Aa = 
1 for different interactions strengths of the BB component 
Ab = 0; 1/2; 1;2; 10 as a function of momentum in units of 
Fermi momentum k/kp. Compare the data to the inset of 

Fig.m 



is greatly increased when the other component forms a 
quasicrystal (see Ab = 10 case in Fig. [8]). 



C. A mixture of Calogero-Sutherland gases with 
interspecies interaction: off-diagonal properties 
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FIG. 9: Inset: projected two-body density matrix P2Aa(-R)i 
Eq. dill) for iVA = A'^B = 10. Aa = Ab = 
0; 0.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9; 1 (from upper to lower 
curves). Main figure: logarithmic scale. Inset: close up of the 
tail. 



the cut of the two-body density matrix for the displace- 
ment R = \xi — x'll = \x2 — x'2\ = L/A as a function 
of the dimer size \x2 — xi\ in Fig. |31 For small values 
of A the function is significantly different from zero and 
the decay is exponentially fast. According to the spec- 
tral decomposition (|15p the function shown in Fig. [3] is 
proportional to the the square of the molecular orbital. 
Of course, this relation has no meaning when ODLRO is 
absent and p2{r,L/2) can become negative for large A. 



Examples of the typical behavior of the one-body den- 
sity matrix in a balanced system tia = are shown in 
Fig.[l]for different values of interaction parameter A. One 
sees how the fermionic nature of AB interactions mani- 
fests in oscillating behavior of pi{x — x') — {'if^x)'i'{x')) . 
While removing a particle from x and moving it to x', 
the ground state wave function changes its sign each time 
x' crosses a particle of the other component. In the non- 
interacting case, A = 0, there are no correlations between 
particle positions and for two particles there is a linear de- 
cay Pia{x) — 7ia(1 — 2nB|a;|/A^B)- In the thermodynamic 
limit the decay in piA(a^) is exponentially fast according 
to Eq. p3|) . but still the one-particle density matrix re- 
mains positive. Instead, for a finite A the interparticle 
correlations are nontrivial and the sign alternations due 
to FTG interactions lead to regions of negative sign. 

We study the ODLRO in the two-body density matrix. 
Fig. ini shows its behavior on the displacement distance R 
in a system with Na = ^Vb = 10 particles and different 
values of interaction parameter. For A = the function 
goes to a constant value as \R\ oo {p2{R) — > 1/20 in 
the present case). The ODLRO gradually vanishes as the 
interaction parameter is increased and at a certain point 
the oscillatory behavior dominates (see inset of Fig. [S]) . 
In the regime where the ODLRO is still present we study 
the shape of the dimer wave function (orbital) . We show 



D. A mixture of Calogero-Sutherland gases with 
interspecies interaction: phase diagram 

In order to understand if the system is superfluid we 
note that the mapping applies also to the excited states. 
The excitation spectrum in a single component Calogero- 
Sutherland model for small frequencies lo can be analyzed 
from the Luttinger liquid theory. The dynamic form fac- 
tor cr(w, k) in the point where the excitation spectrum 
touches zero, k = 2?™, has a power-law dependence on 
the frequency (j{u}, 27m) cx uj^/^~'^ for small w [l^ll^. In 
the weakly interacting regime (i.e. small A) the weight is 
vanishingly small and the behavior of the system is anal- 
ogous a three-dimensional superfluid. Instead, as A is 
increased the system starts behaving as a normal system 
and an impurity drag ged through the system will cause 
an energy dissipation [27[. 

We summarize the phase diagram of a system with 
Aa — Ab in Fig. 1101 The contour plot shows the maximal 
value of the momentum distribution nA{k) as a function 
of A and density imbalance n^/nA- In the noninter act- 
ing limit A — >■ we recover the phase diagram shown in 
Fig. 21 The maximum observed for small density of the 
other component tt-b / tla ^ 1 corresponds to a very high 
peak of the momentum distribution at fc = 0, which we 
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FIG. 10: Altitude of the highest peak in the momentum dis- 
tribution maxn(fc) for A'a = 10 and A^b = 2 — 100 particles 

k 

shown as a contour plot in variables of the density imbalance 
ub/tia and interaction parameter A. The dashed lines sepa- 
rate different phases and correspond to parameters for which 
the height of the peaks at fc = and fe = fcp is equal. Some 
important features of the phase diagram are observed in this 
plot. True Bose-Einstein condensation (constant ODLRO in 
pi (a;)) of A (B) particles happens for A = and ub ~ 
(ua ~ 0). Small, but finite values of A and ne (?ia) corre- 
spond to quasicondensate of A (B) particles with the highest 
peak positioned at A; « 0. By increasing A the position of the 
highest peak shifts to k — kp which corresponds to the phase 
of Fermi-momentum condensation. The strongest divergence 
at the Fermi momentum is observed for A = 1. By increas- 
ing the value of A further the position of the highest peak 
shifts once more to fc ~ which marks the "classical" regime. 
We stop the dashed line for small ne when kp is comparable 
to the momentum quantization 2n/L. In the ideal gas limit 
Ub ^ the Fermi momentum vanishes as fcp — >■ 0. 



associate with a quasicondensate of A particles. At the 
same time in this regime there is a true condensation of 
AA pairs (molecules). The quasicondensate regime ex- 
tends as well to finite values of A where it is limited by a 
power-law decay in terms of the A variable and by an ex- 
ponential decay in terms of ub/tia- A further increase in 
A leads to stronger interparticle correlations and around 
A = 1 we find an appearance of another maximum, which 
this time is situated at the Fermi momentum k — kp. We 
associate this regime with the Fermi-momentum quasi- 
condensation. In order to estimate a border line between 
the different phases we plot with black circles the points 
on the phase diagram where the height of the peak at 
fc = is equal to the height of the peak at fc = fcp- 
For even stronger correlations, A ^ 1, the position of 
the peak shifts again to A; = (as an example momen- 
tum distribution for riA — in Figs. [S]and[S]). In this 
regime the interactions quickly destroy the coherence (see 



A = 10 case in Fig. [T]) and from this point of view the 
systems behaves quite similarly to a classical one. The ki- 
netic energy is much smaller than the repulsion between 
the particles and this leads to creation of a crystal-like 
structure. The momentum distribution can be closely de- 
scribed by a Gaussian with the peak at fc = 0, Eq. (|18p . 
We refer to this regime as a "classical" phase. 

It is important to note that there are no phase transi- 
tions in our system. Instead the dashed lines in Fig. [TU] 
separate different physical regimes and the transition be- 
tween them is continuous (crossover). Furthermore, the 
position of the lines depends explicitly on the number of 
particles. This can be understood by noticing that di- 
vergence at fc = fcp in a Fermi quasicondensate is weaker 
compared to the divergence at fc = in a quasiconden- 
sate; see Fig. [T] 

As discussed above, the system is normal for the value 
A = and arbitrary densities n^^/n-Q, which is an ex- 
cluded line. Instead, for finite but small values of A the 
system behaves practically as a superfluid, while at larger 
values of A there is a strong friction in the system and 
its behavior is normal. In the limit of large A the in- 
terparticle correlations are extremely strong, which is a 
one-dimensional analogue of a crystal. In this regime the 
system is normal and no Bose condensation of any kind 
is present. 



IX. CONCLUSIONS 

To summarize, we have studied properties of a bosonic 
two-component one-dimensional system at zero temper- 
ature. The interactions between particles of the same 
component (AA or BB) are taken to be of a Calogero- 
Sutherland type A(A — '^)/(xi — XjY . The interaction pa- 
rameters Aa and Ab can be changed independently and 
are related in a simple way to the Luttinger parameters 
Ka = 1/Aa, K-q = 1/Ab. In this way the long-range 
properties of the considered model are universal for Lut- 
tinger liquids. For the AB interaction we consider p-wave 
scattering channel with attractive interaction such that 
the wave function of a two body zero-energy scattering 
solution is fiat outside of some short range, i.e., the usual 
fermionic Tonks- Girardeau (FTG) potential. We write 
explicitly the many-body wave function which we obtain 
by a mapping procedure. The same mapping applies to 
the excited states of the system. The phase diagram of 
the system is very rich and is governed by the interac- 
tion strengths Aa, Ab and the density ratio Na/N-q- The 
limiting case of Aa = Ab = corresponds to a mixture 
of two ideal Bose gases with FTG interaction between 
them and was recently studied in Ref. [11]. In this case a 
single component A'b — ideal gas creates a completely 
"true" Bose-Einstein condensate with all particles being 
in the same state. Presence of the other component de- 
stroys the off-diagonal long-range order in the one-body 
density matrix. At the same time there is a divergent 
peak in the momentum distribution at fc = 0, which we 
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interpret as a quasi-condensation. The one-body density 
matrix, two-body density matrix and the momentum dis- 
tribution of such a system can be obtained analyticahy. 
There is a full condensation of AA pairs (molecules) as 
manifested by the presence of the off-diagonal long range 
order in the two-body density matrix, which for large dis- 
placement of a pair factorizes to the product of molecular 
orbitals. There is a BEC-BCS-like crossover from AA- 
pair BEC when ns ^ «a, implying tightly bound AA 
pairs, to extended and strongly overlapping A A Cooper 
pairs when <^ nj^. 

The properties of the system with finite values of A 
are obtained numerically by sampling the known ground 
state wave function with the Monte Carlo technique. The 
interactions between the particles of the same species de- 
stroy the ODLRO both in one- and two- body density 
matrices. Still, the regime of quasicondensation of A par- 
ticles persists for small interaction strength A and small 
concentration of the other component n-Q/nx- The p- 
wave interactions leads to appearance of fermionic-like 
properties, such as oscillations around zero in the one- 
body density matrix, which leads to a divergence in the 
momentum distribution nx{k) at the Fermi momentum 
of the other component k = ±A;p = ±7rnB. We refer to 
this regime as a phase of Fermi-momentum quasiconden- 
sation. 

By making the interactions with the species stronger 



we reach the regime where the potential energy of repul- 
sive interaction is much larger than the kinetic energy 
which leads to appearance of a crystal-like order. We re- 
fer to this regime as "classical". There is a very rapid 
decay of coherence in the one-body density matrix. The 
momentum distribution is well approximated by a Gaus- 
sian, for which we provide an explicit expression. 

Finally, superfluid properties of the system are dis- 
cussed on the basis of the Luttinger liquid approach. 
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